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ABSTRACT 

Aims. We study the magnetic field structure, strength, and energy density in the Scd galaxy M33. 
Methods. Using the linearly polarized intensity and polarization angle data at 3.6, 6.2 and 20 cm, we determine 
variations of Faraday rotation and depolarization across M33. We fit a 3-D model of the regular magnetic field to 
the observed azimuthal distribution of polarization angles. We also analyze the spatial variation of depolarization 
across the galaxy. 

Results. Faraday rotation, measured between 3.6 and 6.2 cm at an angular resolution of 3' (0.7 kpc), shows more 
variation in the south than in the north of the galaxy. About 10% of the nonthermal emission from M33 at 3.6 cm 
is polarized. High degrees of polarization of the synchrotron emission (> 20%) and strong regular magnetic fields 
in the sky plane (~ 6.6 fiG) exist in-between two northern spiral arms. We estimate the average total and regular 
magnetic field strengths in M33 as ~ 6.4 and 2.5 /iG, respectively. Under the assumption that the disk of M33 is 
flat, the regular magnetic field consists of horizontal and vertical components: however the inferred vertical field 
may be partly due to a galactic warp. The horizontal field is represented by an axisymmetric (m = 0) mode from 
1 to 3 kpc radius and a superposition of axisymmetric and bisymmetric (m = + 1) modes from 3 to 5 kpc radius. 
Conclusions. An excess of differential Faraday rotation in the southern half together with strong Faraday dispersion 
in the southern spiral arms seem to be responsible for the north-south asymmetry in the observed wavelength 
dependent depolarization. The presence of an axisymmetric m — mode of the regular magnetic field in each 
ring suggests that a galactic dynamo is operating in M33. The pitch angles of the spiral regular magnetic field 
are generally smaller than the pitch angles of the optical spiral arms but are twice as big as simple estimates 
based on the mean-field dynamo theory and M33's rotation curve. Generation of interstellar magnetic fields from 
turbulent gas motions in M33 is indicated by the equipartition of turbulent and magnetic energy densities. 

Key words, galaxies: individual: M33 - radio continuum: galaxies - galaxies: magnetic field - galaxies: ISM 



1. Introduction 

Magnetic fields in galaxies can be traced by radio 
polarization measurements. The linearly polarized 
intensity gives information about the magnetic field 
strength and orientation in the plane of the sky, 
Faraday rotation measurements enable us to de- 
termine the magnetic field strength and direction 
along the line of sight and depolarizing effects can 
be sensitive to both components. 

M33, the nearest Scd galaxy at a distance of 
840 kpc, with its large angular size (high-frequency 
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radio contin uum emission was detect ed in an area 
of 35' X 40', iTabatabaei et al.ll2007cl ) and medium 



inclination of 56°, allows determination of the mag- 
netic field components both parallel and perpen- 
dicular to the line of sight equally well. The first 
RM study of M33, based o n polarizati on observa- 
tions at 11.1 and 21.1cm suggested 
a bisymmetric r egular magnetic field s tructu re in 
the disk of M33. iBuczilowski fc Beck I (|l99ll ) con- 
firmed the presence of this bisymmetric field using 
two further polarization maps at 6.3 and 17.4 cm. 
However, these results may be affected by weak 
polarized intensity and the consequent high uncer- 
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tainty in RM in the southern half of M33, due to the 
low-resolution (7.7' or 1.8 kpc) and low-sensitivity 
of the observations. 

Our recent observations of this galaxy pro- 
vided high-sensitivity and high-resolution (3' ~ 
0.7 kpc) maps of total power and linearly polar- 
ized intensity at 3.6, 6.2 , and 20 cm presented by 



Table 1. Integrated nonthermal flux densities and 
average nonthermal degree of polarization at 180". 



( Tabatabaei et al.ll2007cf ). These data are ideal to 



study the rotation measure (RM), magnetic fields 
(structure and strength) , and depolarization effects 
in detail. 



( Tabatabaei et al.l l2007cr ) found a north-south 
asymmetry in the polarization distribution that 
is wavelength-dependent, indicating a possible 
north-south asymmetry in Faraday depolarization. 
Investigation of this possibility requires a knowl- 
edge of the distribution of RM, turbulent magnetic 
field and thermal electron density in the galaxy. 
Furthermore, depolarization is best quantified using 
the nonthermal degree of polarization rather than 
the fraction of the total radio em ission that is po- 
larized. Tabatabaei et al.l (j2007br ) developed a new 
method to separate the thermal and nonthermal 
components of the radio continuum emission from 
M33 which yielded maps of the nonthermal syn- 
chrotron emission and the synchr otron spectral in 



dex va riations across the galaxy ([Tabatabaei et al 



2007al) . The nonthermal maps are used in this paper 
to determine the nonthermal degree of polarization 
at different wavelengths, and hence to study how 
the radio continuum emission from different parts 
of M33 is depolarized. 

By combining an analysis of multi-wavelength 
polarization angles with mo deling of the wav e length 
dep endent depolariza t ion, Fletcher et al.l (|2004l ) 



and iBerkhuiisen et al.l (|l997l ) derived the 3-D reg- 



ular magnetic field structures in M31 and M51, re- 
spectively. The high sensitivity of our new observa- 
tions allows a similar study for M33. 

We first determine the nonthermal degree of po- 
larization using the new nonthermal maps (Sect.[2|). 
Then we calculate a map of the RM intrinsic to 
M33 with a spatial resolution of 3' or 0.7 kpc and 
probe its mean value in rings in the galactic plane 
in Sect. [3l The regular magnetic field structure is 
derived in Sect. |4] and the estimated strengths are 
presented. We derive a map for the observed depo- 
larization and discuss the possible physical causes of 
depolarization sources in Sect. [51 Furthermore, we 
discuss the resulting vertical fields and pitch angles 
in Sect. 6. The estimated energy density of the mag- 
netic field is compared to the thermal and turbulent 
energy densities of the interstellar gas. 



A 


Snth 


Spi 


Pnth 


(cm) 


(mJy) 


(mJy) 


% 


3.6 


370 ± 60 


38±4 


10.3 ±2.0 


6.2 


696 ±110 


79 ±5 


11.3± 1.9 


20 


1740 ± 65 


115 ± 10 


6.6 ±0.6 



2. Nonthermal degree of polarization 

The degree to which synchrotron emission is polar- 
ized reflects the degree of coherent structure in the 
magnetic field within one beam purely reg- 

ular magnetic field will produce about 759([J frac- 
tional polarization of synchrotron emission. The 
quantity of interest is the degree of polarization of 
the synchrotron emission or 'nonthermal degree of 
polarization', Pnth = PI/-^nth, where PI is the in- 
tensity of the linearly polarized emission and /nth 
is the intensity of the nonthermal emission. 

Since we observe the total intensity /, which 
consists of both nonthermal /nth cind thermal /th 
intensities at cm wavelengths, Pnth cannot be cal- 
culated straightforwardly. To date, Pnth has been 
estimated by assuming either a fixed ratio Ith/I 
or /nth has been derived assuming a fixed spectral 
index of the synchrotron emission. Here we use a 
new, more robust method for determining the dis- 
tribution of /nth by correcting Ha maps for dust 
extinction, using multi-frequency infra-red maps at 
the same r esolution, and thus indep endently esti- 
mating Lh ([Tabatabaei et al.l[2007bi) . Using the PI 
maps of [Tabatabaei et all tOOld) and the nonther- 
mal maps obtained bv [Tabatabaei et"aL (|2007bh . 
we derived maps of Pnth at different wavelengths. 

Figure [1] shows Pnth at 3.6 cm. High non- 
thermal degrees of polarization (Pnth > 30%) 
are found in several patches of M33, with the 
most extended region of high Pnth in the north- 
ern part of the magne tic filament identified by 
Tabatabaei et al. I ([2007cl) . inside the second contour 
at DEC > 30° 54' in Fig.[Il The high Pnth also exist 
at 6.2 cm in these regions. 

Integrating the polarized and nonthermal inten- 
sity maps in the galactic plane out to galactocentric 
radius of R < 7.5 kpc, we obtained the flux densities 
of the nonthermal Snth and linearly polarized Spi 
emission along with the average nonthermal degree 
of polarization Pnth = Spi/Snth- Table [Tj gives Snth, 
Spi, and Pnth at different wavelengths, all at the 
same angular resolution of 180". At 3.6 and 6.2 cm, 
Pnth is the same, demonstrating that Faraday depo- 



^ If a nonthermal spectral index of a„ = 1 is used. 
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Fig. 1. Nonthermal degree of polarization at 
3.6 cm, with an angular resolution of 2' (the beam 
area is shown in the left-hand corner). Overlayed 
are contours of the linearly polarized intensity at 
3.6 cm with levels of 0.1 and 0.3 mJy/beam. 



larization effects are not significant at these wave- 
lengths: the weaker Spi at 3.6 cm is due to lower 
synchrotron emissivity, as expected from the power- 
law behavior of synchrotron emission with respect 
to frequency. However, Faraday depolarization ef- 
fects are possibly important at 20 cm reducing Pnth- 



3. Rotation measures 

When linearly polarized radio waves propagate in 
a magneto-ionic medium, their polarization vector 
is systematically rotated. The amount of Faraday 
rotation depends on the wavelength (A) of the radio 
emission, the strength of the magnetic field along 
the line of sight (that is the regular field in the line 
of sight Bj-cgii), and the number density of thermal 
electrons (rig) along the line of sight (/): 



rad 



0.81 fA 

\m 



Br, 



1-3 



gle (0=1 arctan^) differs from the intrinsic polar- 
ization angle {(jfi) as 



(1 



where L is the path length through the magneto- 
ionic medium. Hence, the measured polarization an- 



^ A<j) = \^ n. 



(2) 



When Faraday depolarization is small (Faraday- 
thin condition), TZ d oes not depend on wavelength 
(ISokoloff et al.lll998h and an estimate for TZ can be 
obtained from measurements of the polarization an- 
gles at two wavelengths: 



RM 



radi 



(01 /rad) - ((/'2/rad) 



(3) 



In this definition, the unknown intrinsic polariza- 
tion angle of the source (or sources along the line 
of sight) cancels. Positive RM indicates that -Brcg|| 
points towards us. 

Part of the measured RM is due to the interstel- 
lar medium of M33 (intrinsic RM or RMj), the rest 
is due to the Galactic foreground medium (RMfg), 
RM = RMi + RMfg. The foreground rotation 
measure in the direction of M33 is mainly caused 
by the extended Galactic magnetic bubble identi- 
fied as region A bv lSimard-Normandin fc Kronberg 
Assuming that the intrinsic contribu- 
tions of the extragalactic sources 3C41, 3C42, 
and 3C48 (in a 5° x 5° region around M33) 
themselves cancel out and the intra g alactic con- 
tribution is neg l igible , Broten et al. ( 19881 ) and 
Tabara fc Inouel (|l980l ) found a foreground rota- 
tion measure of — 57 ± 10 rad m~^ for those sources. 
For the polarized sou rc es in the 2° x 2° M33 field. 



Buczilowski fc Beck I (|199l[ ) found a foreground 
RM of —55 ± 10 rad m~^. About the same v alue 
was obtained bv iJohnston-Hollitt et al.l (|2004[) . In 
the following we use RMfg = — 55radm~^. 



Using the polarization data of iTabatabaei et al 



(2007g), we first obtained the distribution of RM 
between 3.6 and 20 cm across M33 (Fig. [21 left 
panel), showing a smooth distribution of RM in 
the northern half of the galaxy. However, stronger 
and more abrupt fluctuations in RM occur in the 
southern half of the galaxy, which are not due to 
the ±n73radm~^ ambiguity in RM between these 
wavelengths (±n7r/ | — |). Weak polarized 
em ission in the southern h alf at 20 cm (presented 



Tabatabaei et al.ll2007cl ) can be finked to these 



RM variations. Between 3.6 and 6.2 cm, RM varies 
less than between 3.6 and 20 cm in the south of 
the galaxy (Fig. [21 right panel). This indicates that 
the relation between A0 and A~^ in Eq. ^ is not 
linear over the interval between 3.6 and 20 cm due 
to Faraday depolarization at 20 cm in the south of 
M33 and so RM measured at these wavelengths is 
not a good estimator for TZ. 
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Fig. 2. Left: observed rotation measure map of M33 (radm~^) between 3.6 and 20 cm with contours of 
3.6 cm polarized intensity. Contour levels are 0.1, 0.2, 0.4, 0.6, 0.8mJy/beam. Right: observed rotation 
measure between 3.6 and 6.2 cm with contours of 6.2 cm polarized intensity. Contour levels are 0.3, 0.4, 
0.8, 1.2, 1.6 mjy/beam. The angular resolution in both maps is 3' (the beam area is shown in the left-hand 
corners). The straight line in the left panel shows the minor axis of M33. 



Figure [3] shows RMi between 3.6 and 6.2 cm 
which varies in a range including both positive and 
negative values. Comparing RMj with the overlayed 
contours of P^, RMj seems to vary more smoothly 
in regions of high PI. The apparent agreement be- 
tween the ordered magnetic field in the plane of the 
sky and the regular magnetic field in the line of sight 
is clearest in the north and along the minor axis and 
is best visible in the magnetic filament between the 
arms IV and V in the north-we st of M33 (Fig. 01 
see also Tabatabaei et al.|[2007cl) where RMi shows 
small variation within the PI contours. This indi- 
cates that the ordered magnetic field in this region 
is mainly regular. Sign variations of RM.^ are more 
frequent (arising locally) in the southern half, where 

^ Note that PI is related to the magnetic field in the 
plane of the sky that is a combination of both a mean 
field (or coherent regular field) and anisotropic (com- 
pressed or sheared) random fields. RMi is related to 
only coherent regular field along the line of sight. 



there is no correlation with PI, than in the northern 
half of the galaxy. This indicates that the regular 
magnetic field is more affected by local phenomena 
(like s tarforming activity, e.g. see iTabatabaei et al] 
l2007d) in the south than in the north of the galaxy. 
The local RM^ variations between large positive 
and negative values may represent loop-like mag- 
netic field structures going up from and down to 
the plane (e.g. Parker loops). This is particularly 
seen in the central part of the galaxy besides re- 
gions in the southern arm IIS (Fig. [5]). 

Figure [3] also shows that the magnetic field is 
directed towards us on the western minor axis (at 
azimuth 6* R:! 110° and 9 w 290°), but has an op- 
posite direction on the eastern side. The large RMj 
values in regions with small electron density, e.g. on 
the eastern and western minor axis and in a clumpy 
distribution in the central sout h near the major 
axis with 30° 25' <DEC<30° 30' (ITabatabaei et all 



2007br ). indicates a strong magnetic field along the 
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Fig. 3. Left: rotation measure map of M33 (radm"'^) between 3.6 and 6.2 cm after correction for the 
foreground RMfg = — 55radm~^, with an angular resolution of 3' (the beam area is shown in the left- 
hand corner). Overlayed are contours of 6.2 cm polarized intensity. Contour levels are 0.3, 0.4, 0.8, 1.2, 1.6, 
3.2, 6.4mJy/beam. Also shown are rings at 1, 3, and 5kpc radii as used in Sect.[?l Right: the distribution 
of the estimated error in RM is shown in the right panel. 



line of sight and/or large path length through the 
magneto-ionic medium. Considerable Faraday ro- 
tation measured on the eastern and western mi- 
nor axis hints to deviation from a purely toroidal 
struc t ure for the l arge-s cale magnetic field (iKrause 
19901 : iBeck et al.l[l996l) . Furthermore, along with 
the magnetic field parallel to the disk, the presence 
of a vertical field component to the galactic disk 
is indicated in kpc-scale regions of large RM^ but 
small Ue values. The existence of the vertical mag- 
netic field in M33 (that is strong near the major 
axis) is shown in Sect. HI 



4. The magnetic field 

4.1. The regular magnetic field structure 

The A3. 6 cm polarization angles were corrected for 
Faraday rotation and rotated by 90° to obtain the 
intrinsic orientation of the regular magnetic field 
component in the plane of the sky ( Breg±)- Figure[3] 
shows the derived i?rcg_L superimposed on an opti- 
cal image of M33. The orientation of Bycg± shows 
a spiral magnetic field pattern with a similar pitch 
angle to the optical arms in the north and south. 



but larger pitch angles in the east and west of the 
galaxy. 

The apparent continuity of the regular magnetic 
field straight through the center of M33 is remark- 
able. The Faraday rotation map (Fig. [3]) shows that 
the field is oppositely directed on the east and west 
sides of the center. Unfortunately the current data 
lack the resolution to investigate the field properties 
in this region further. 

The behavior of -Brcg± along the minor axis 
{9 « 110°, e w 290°) is also informative. B^g^ 
on the minor axis clearly has a strong radial com- 
ponent (Fig. [4]). If we assume that i?rcg lies solely in 
the disc, so that -Brog|| is produced by the galaxy's 
inclination (i.e. = 0), then the change of sign 
in RM along the minor axis (Fig. [3]) indicates that 
the direction of the radial component of Breg is to- 
wards center on both sides of the minor axis. This 
means that the dominant azimuthal mode of i?rog 
cannot be the bis ymmetric 171 = 1 mode suggested 
by ea rlier studies ( Beck|[l979t Buczilowski fc Beck 



I l99l| ): we expect the dominant mode to be even. 
In order to identify the 3-D structure of the 



regular magnetic field 



fit a parameter- 
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Fig. 4. Left: Optical image (B-band, taken from the STScI Digitized Sky Survey) of M33 with overlayed 
vectors of intrinsic magnetic field in the sky-plane with 3' angular resolution. The length gives the 
polarized intensity at 3.6 cm, w here 1' is equivalent to 0.37mJy/beam. Right: A sketch of the optical arms 
(jSandage &: HumphrevjilOSOl) overlayed on the linearly polarized in tensity with an angular resolution of 
2' (contours and grey scale) at 3.6 cm (see Tabatabaei et ahlboOTcl) . 



ized model of -Brcg to the observed polarization 
angles at different wavel engths using the method 
succ essfully employ e d by iBerkhuiisen et al.l (|l99l1 ) 
and iFletcher et al.l (|2004l ) to determine the mag- 
netic field structures of M51 and M31 respectively. 



At each wavelength the maps in the Stokes pa- 
rameters Q and U are averaged in sectors of 10° 
opening angle and 2 kpc radial width in the range 
1 < r < 5 kpc. The size of the sectors is chosen 
so that the area of the smallest sector is roughly 
equivalent to one beam-area. We also take care 
that the standard deviation of the polarization an- 
gle in each sector is greater than the noise (if the 
sector sizes are too small then fluctuations in an- 
gle due to noise can be larger than the standard 
deviation). Then the average Q and U intensities 
are combined to give the average polarization an- 
gle i/j — 0.5 arctanC//Q and the average polarized 
emission intensity Pl^/{CP + lP) in each sector. 



The observed polarization angles are related to 
the underlying properties of the regular magnetic 
field in M33 bji 



^ = MB.eg±) + A2RMi(B,eg||) + A^RMfj 



(4) 



where ipo is the intrinsic angle of polarized emission, 
RMi is the Faraday rotation experienced by a pho- 
ton as it passes through the magneto-ionic medium 
of M33 and RMfg is foreground Faraday rotation 
due to the Milky Way. i?rcg± is the component of 
the regular magnetic field of M33 that lies in the 
sky-plane and -Brcg|| the component directed along 
the line-of-sight. 

The cylindrical components of the regular field 
in the disk of M33 .Breg = {Br, Be, Bz) can be 
represented in terms of the Fourier series in the az- 
imuthal angle 9: 

Br = Bosinpo + BismpiCos{9 ~ (3i) 

Note that we choose different symbols for the aver- 
age polarization angle <f) and the average polarization 
angle in sectors tp. 
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+ 3-2 sinp2 cos 2(6' - /32), 

i?o cospo + -Bi cospi cos(0 — /3i) (5) 

+ B2 C0SP2 cos 2(6* - /32), 

B.o + B,i cos{0 - (3,i) + B,2 cos 2(0 - /3,2), 



where Bm and Bzm are the ampUtude of the mode 
with azimuthal wave number m in the horizontal 
and vertical fields, Pm is the constant pitch angle 
of the m'th horizontal Fourier mode (i.e. the an- 
gle between the field and the local circumference) 
and Pm and j3zm are the azimuths where the non- 
axisymmetric modes are maximum. The amplitudes 
of the Fourier modes are obtained in terms of the 
variables Bm in units of rad m~^ : in order to obtain 
amplitudes in Gauss independent information is re- 
quired about the average thermal electron density 
and disc scale-height. Useful equations describing 
how -Brog± and i3rog|| are related to the field com- 
ponents of Eg. ^ can b e found in Appendix A of 



Berkhuiisen et al.1 (|1997[ ). 

Using Eqs. @ and (O, and taking into ac- 
count the inclination (56°) and major-axis orien- 
tation (23°) of M33, we fit the three-dimensional 
-Breg to all of the observed polarization angles in a 
ring by minimizing the residual 



(6) 



where ^„ is the observed angle of polarization, 
tp{On) the modelled angle in the sector centred on 
azimuth 0„ and fT„ are the observational errors. 
The Fisher test is used to ensure that the fits at 
each wavelengt h are statistically, equall y good (see 
Appendix B in Berkhuiisen et al]|l997[) . Errors in 
the fit parameters are estimated from the region of 
parameter space where 5* < ^t the 2a level. 

4.2. Results of fitting 

We applied the method described in Sect. 2] to the 
polarization maps at 3.6 and 6.2 cm. A preliminary 
examination of the data at all three wavelengths 
showed that the 20 cm polarization angles do not 
have the same oa Faraday rotation depen- 
dence as the angles at 3.6 and 6.2 cm. This is most 
probably because the 20 cm signal is strongly de- 
polarized by Faraday effects (see Sect. [5]) and so 
only photons from an upper layer of the emitting 
region are detected in polarization. If the depo- 
larization is constant at a given radius, a 'depth' 
para meter can be used to take account of this ef - 



fect (jBerkhuiisen et al.lll997t iFletcher et al.ll2004) . 



However, in the case of M33 depolarization is 
strongly asymmetric (Fig. llip and this method does 




(a) 
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Fig. 5. Polarized intensity (PI) at 6 cm in the ring 
1-3 kpc, observed (points with error bars) and ex- 
pected from the modelled regular magnetic field 
(lines) with PIcx B^^^j_ given by the different fits 
shown in Table. [H (a) using a fitted regular mag- 
netic field containing the modes m = + zO + zl. 
(&) as (a) but for the fit containing the modes 
m = -f 1 -I- zl. (c) as (a) but for the fit containing 
the modes to = + 2. 



not lead to consistent results: we therefore work 
only with the polarization angles at 3.6 and 6.2 cm 
to model the regular magnetic field. 

We fixed the foreground Faraday rotation, 
RMfg in Eq. (gl), to -SSradm"^ (see Sect.E]). For 
each ring we found more than one statistically good 
fit to the observed polarization angles, three differ- 
ent fits in the ring 1-3 kpc and two fits in the ring 
3-5 kpc, the fitted parameters are shown in TableO 
All of the fits require the presence of more than 
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50 100 150 200 250 300 350 
azimuth [degrees] 

Fig. 6. As figure [5] but for the ring 3-5 kpc using a 
fitted regular magnetic field containing the modes 
m = 0+l + zl (top) and m = + 2 {bottom). 



one azimuthal Fourier mode and have two common 
characteristics: the presence of an m = mode that 
has a significant amplitude; the pitch angle of the 
m — mode is in the range 40° ^ p ^ 50°. The 
reason why several equally good fits are found is 
the rather weak large-scale intrinsic rotation mea- 
sure signal, i.e. the amplitude of systematic rotation 
measure variations is rather low compared to local 
fluctuations. This is a sign that the regular mag- 
netic field of M33 is not as well-ordered and strong, 
relative to t he small-scale field, a s that of, for ex- 
ample, M31 (jFletcher et al.ll2004h . 

The model regular magnetic field given by Eq. [5] 
is fitted to the observed polarization angles^ in or- 
der to obtain the results shown in Table [2l Since 
we have not made use of the observed polarized 
intensity we can try to use this to select the best 
regular magnetic field model for each of the two 
rings from the fits given in Table [51 We compare 
the predicted azimuthal pattern of polarized inten- 
sity from the model fields with the observed po- 
larized intensity (PI) at 6.2 (Faraday depolariza- 
tion effects are negligible at 6.2 cm and emissivity 
is higher than at 3.6 cm due to the spectral index, 
thus giving a stronger signal) . The model described 
by Eq. [5] is not designed to reproduce the observed 
PI so we cannot make meaningful statistical assess- 
ments about the relative merits of the different fits 



in a given ring. But we can judge whether or not the 
fits are better or worse than each other in explain- 
ing the location of the main maxima and minima 
in the observed pattern of PI and so try to select a 
preferred model field for each ring. 

Figures [5] and [H] shows the square of the pre- 
dicted plane-of-sky regular magnetic field S^egj. 
for each of the fits given in Table [2] and the ob- 
served 6.2 cm polarized intensity, both normalized 
to avoid having to use a prescription for the poorly 
known synchrotron emissivity. In the case of en- 
ergy equipartition between cosmic rays and mag- 
netic fields the polarized intensity would be pro- 
portional to a higher power of i3rcg± than 2 and 
maxima would be more pronounced. 

In the ring 1-3 kpc, the model field with the 
components ofm = 0-|-zO-|-zlQ (Fig. [5^) repro- 
duces the broad features of the observed polarized 
emission better than the models using the modes 
■m = 0+l + zl (Fig. [5b) and m = -I- 2 (Fig. [5fc). 
The match to the observed PI is far from perfect 
in Fig. [5^ but this is the only model field that can 
account for the strong excess of PI at 6 cm in the 
northern half of the disc at these radii. In the ring 
3-5 kpc the fit using m = -I- 1 -I- zl (Fig. [H top) 
is better at reproducing the general pattern of po- 
larized emission at 6 cm than the other statistically 
good fit using m = 0-1-2 (Fig.[6l bottom). Again, the 
match to observations is not perfect, but the model 
with m = + 2 would produce a strong maximum 
in PI at 6* ^ 90° that is not observed. 

To summarize: we select the statistically good fit 
using the modes m = + zO + zl in the ring 1-3 kpc 
and that using m — 0-|-l-|-zlin the ring 3-5 kpc as 
being the best descriptions of the regular magnetic 
field in M33 (the parameters of these two preferred 
fits are given in columns 3 and 6 of Table [2]). Our 
reason is that the B'^egj_ produced by these models 
produces a much closer match to the observed pat- 
tern of PI at 6.2 cm than other statistically good 
models. Fig. [7| shows the regular magnetic field in 
a face-on view of the galaxy (thus the vertical field 
components cannot be seen). 

4.3. The equipartition magnetic field strengths 

The strengths of the total magnetic field Btot and 
its regular component i?rog can be found from the 
total synchrotron intensity and its degree of linear 
polarization Pnth- Assuming equipartition between 
the energy densities of the magnetic field and cos- 
mic rays {ecR = estot = ^tot/87r), 

* zO and zl are the first and secound Fourier modes 
of the vertical field. 
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Table 2. Parameters of the fitted models and their 2cr errors. RM[g is the Faraday rotation measure 
arising in the Milky Way, and Pm are the amplitude and pitch angle of the mode with wave number 
m, and /3„j is the azimuth where a mode with azimuthal wave number m is maximum. The minimum 
value of the residual S and the value of are shown for each model in the bottom lines. The combination 
of azimuthal modes of m = + zO + zl and m = + 1 + zl can best re-produce the observed polarized 
intensity in the 1-3 and 3-5 kpc rings, respectively (Figs and [S]) . 



Units Radial range (kpc) 
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O + zO + ^l 
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+ 1 + zl 


+ 1 + zl 
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56 


83 


57 


57 


82 


x' 






85 


85 


84 


84 


85 



(2a„ - l)c2{a„)Lc3 



^(7) 



(|Beck fc Krausell2005l ). where K' = K + 1 with K 
the ratio between the number densities of cosmic 
ray protons and electrons, /„ is the nonthermal in- 
tensity, L the pathlcngth through the synchrotron 
emitting medium, and an the mean synchrotron 



spectral index. Ep = 938.28MeV: 
is the proton rest energy and 



Cl 

C2(an) 



3e/{4TTmJe') 
1 



6.26428 • IQis 



erg^.s. G 

(a„ + 5/3) /{an + 1) r[(3a„ + l)/6] 
xr[(3a„ + 5)/6]. 



For a region where the field is completely regu- 
lar and has a constant inclination i with respect 
to the sky plane {i = 0° is the face-on view), 
C3 = [cos (1)]*^""+^' . If the field is completely tur- 
bulent and has an isotropic angle distribution in 



three dimensions, C3 — (2/3)^""+^-'/^. If the syn- 
chrotron intensity is averaged over a large volume, 
[cos (i)](""+^) has to be replaced by its average over 
all occurring values of i. 

The strength of the regular magnetic field in the 
plane of the sky can be estimated f rom the observed 
nonth ermal degree of polarization (jSegalovitz et al 



1.50 X 10 ^ erg p 



nth 



37 



37 + 7 



^ (l-g)^V^r[(7 + 5)/4] 
2qr[{j + 7)/4]F{z) 



F{i) 



271 



ri-sin^sin2^)<^+'^/' de, 



(8) 



with Srcg/Btur = g2/(i+7)^ J ^ 2an + 1, and 6 
the azimuthal angle ( Btur is the turbulent mag- 
netic field). This formula assumes that the regular 
magnetic field has a single orientation, is parallel to 
the disk and, taken over the galaxy as a whole, has 
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The equipartition total magnetic field strength 
of 6 ± 0.5^G is s lightly higher than the to- 



Fig. 7. The regular magnetic field described by our 
favoured fitted model (see Sect. l4.2l for details). The 
galaxy has been deprojected into a face-on view, 
so only the disk-plane components of the magnetic 
field are shown, the vertical components are not 
visible. The dashed line shows the major axis of 
M33, with north to the top. The ring boundaries 
are at 1, 3 and 5 kpc. 



no further preferential orientation with respect to 
any fixed direction in space. 

The determined average values of /„, q;„, and 
Pnth with the assumed values of K (~ 100) and L (~ 
1 kpc/cos«) lead to Btot ~ 6.4 ± 0.5 /iG and Bj-og — 
2.5 ± l.O^JLG for the disk of M33 {R <7.5kpc). The 
strongest regular magnetic field is found in between 
the northern arms IV N and VN (in the magnetic 
filament) with Breg — 6.6 /iG where B^ot — 8.3 /xG. 

The regular magnetic field strength estimated 
from the mean rotation measure in sectors and as- 
suming Tie ~ 0.03cm~3 (see Sect. [5|) is 1.4 l^J /iG 
and 0.6 ± 0.1 /xG in the ring 1-3 kpc and 3-5 kpc, 
respectively. In the ring 1-3 kpc, the regular field 
strength is consistent with that estimated from the 
equipartition assumption, while it is much smaller 
in the ring 3-5 kpc. The more frequent RM varia- 
tions in both amplitude and sign in the second ring 
(see Fig. 3) indicate a line-of-sight field component 
with many reversals and a small mean RM within 
each sector, whereas the field component in the sky 
plane causes significant polarized emission (which 
is insensitive to the field reversals). 



tal field obtained by iBuczilowski fc Beck I (|l99lh . 
4 ± 1/iG, assuming the minimum total energy re- 
quirement of the disk-like synchrotron source. This 
is not surprising because a difference of ~20% 
is expected between the minimum and equiparti- 
tion magneti c field strengths for fi eld strengths of 
about 5 /iG (jBeck fc Krausd 120051 ). Furthermore, 
the mean nonthermal fraction from the standard 
thermal/nonthcrmal separation method is lower 



than t hat from the new method ([Tabatabaei et al 
2007bl ). resulting in a weaker equipartition mag- 
netic field. 

As the polarized intensity (PI/0.75, corrected 
for the maximum fractional polarization in a com- 
pletely regular field) is related to the regular mag- 
netic field strength, and the nonthermal intensity 
(/„) to the total magnetic field strength in the 
plane of the sky, /„ — (PI/0.75) gives the nonther- 
mal emission due to the turbulent magnetic field 
-Btur ■ Using this intensity with Eq. ^ and assuming 
a completely turbulent field yields the distribution 
of Btur across the galaxy. Figure [5] shows strong 
-Stur (> 7/iG) in the central region of the galaxy, 
the arm IS, and parts of the northern arm IN. 

Using the mean synchrotron flux density, syn- 
chrotron spectral index, and degree of polarization 
in rings, we also derive the average field strengths in 
rings. Figure [9] shows some fluctuations but no sys- 
tematic increase or decrease of these strengths with 
galactocentric radius. The small bump at 4.5 < i? < 
5.5 kpc is due to the M33's magnetic filament. 



5. Depolarization 

The depolarization observed at a certain wave- 
length is defined as the ratio of the nonthermal de- 
gree of linear polarization Pnth and the theoretical 
maximum value po (75% for Q!„ = 1). Generally, 
depolarization may be caused by instrumental ef- 
fects as the bandwidth and beamwidth of the ob- 
servations or by the wavelength-dependent Faraday 
depolarization. Bandwidth depolarization occurs 
when the polarization angles vary across the fre- 
quency band, reducing the observed amount of po- 
larized emission. It is given by smc(2RMA^ 5v/v), 
where 5v is the bandwidth of the observations 
(e.g. iReichI 120061 ). In our study, the wavelengths. 



bandwidths and RM values lead to a negligible 
bandwidth depolarization. Beamwidth depolariza- 
tion occurs when polarization vectors of different 
orientation are unresolved in the telescope beam. In 
order to compensate this effect, the ratio of the non- 
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Fig. 8. Distribution of tlie turbulent magnetic field 
strength, i?tur (/^G), in M33 with an angular reso- 
lution of 3' (the beam area is shown in the left-hand 
corner). 



thermal degree of polarization at two wavelengths 
is used at a same angular resolution, 



DP 



A2/A1 



nth 
nth 



(9) 



where, A2 > Ai. The observed depolarization 
DP^jj/Aii that is only wavelength dependent, is 
called the Faraday depolarization. 

We derived the depolarization DP20/3.6 using 
the maps of nonthermal and polarized intensity 
at 20 and 3.6 cm at the same angular resolution 
of 180" (Fig. [n top panel). The southern half 
of the galaxy is highly depolarized compared to 
the northern half. While DP20/3.6 changes between 
0.0 and 0.5 in the south, it varies between 0.3 
and 1.0 in the north. Considerable depolarization 
is found at the positions of the prominent HII re- 
gions NGC604 (RA = 1'' 34™ 32. 9^ DEC = 30° 47' 
19.6"), NGC595 (RA = l'' 33™ 32.4", DEC ==30° 
41' 50") and IC133 (RA = l'' 33" 15. 3^ DEC = 30° 
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Fig. 9. Variation of the mean total, regular, and 
turbulent magnetic field strengths in rings of 500 pc 
width with galactocentric radius. 



53' 19.7") as can be expected due to their high den- 
sities of thermal electrons. The strongest depolar- 
ization in the inner galaxy occurs in the main south- 
ern arm I S. No depolarization (DP ~ 1) is seen on 
the eastern end of the minor axis and some northern 
regions. 

There are several mechanisms that can lead 
to wa v elengt h -dependent Fa r aday depolarization 
(|Burnl 119661 : ISokoloff et aP llQQSh . Differential 
Faraday rotation occurs when synchrotron emission 
originates in a magneto-ionic medium containing 
a regular magnetic field. The polarization plane of 
the radiation produced at different depths within 
the source is rotated over different angles by the 
Faraday effect and this results in a decrease in the 
measured degree of polarization. Faraday dispersion 
is depolarization due to fluctuations in the rotation 
measure within a beam, caused by the turbulent 
magnetic field and distribution of thermal electrons 
along the line of sight. When this dispersion is in- 
trinsic to the source, it is called internal Faraday 
dispersion. In case of a dispersion in an external 
screen it is called external Faraday dispersion. This 
depolarization effect may be responsible for the 
north-south asymmetry in the polarized emission 
from M33, if an asymmetry in distribution of the 
foreground magneto-ionic medium exists. However, 
as M33 cannot be resolved in the available fore- 



ground surveys, like RM (jJohnston-Hollitt et al 



20041 ) and Ha (Wisconsin Ha mapper) surveys, we 



do not discuss this depolarization further. Finally 
rotation measure gradients on the scale of the beam 
or larger due to systematic variation in the regular 
magnetic field can also lead to depolarization. The 
regular field in M33 is not strong enough nor is the 
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inclination of the galaxy high enough for this effect 
to be significant (in contrast the highly regular field 
of the strongly incline d galaxy M31 does produce 
strong RM gradients, iFletcher et al.ll20o4 ). 

The depolariza tion due to internal Faraday dis- 
persion (given by ISokolofF et al. is 



CRM 



0.81 (n,) Btur VLd/f, 



(10) 



where the dispersion in rotation measure is (Trm, 
with L the pathlength through the ionized medium, 
/ the filling factor of the Fa raday-rotat ing gas along 



Becljl2007l). and d th e 



Ohno & Shib"atalll993[ ) 



the line of sight (~0.5, 
turbulent scale (~ 50 pc, 
Using the Ha emission measure {EM = J nl.dl) 
and a clumping factor /c = {n^)"^ / {n1) describ- 
ing the variations of the electron density, (rie) can 
be determined by (rie) = fc EM/L. For the 
local interstellar medium, [Manchester fc Meboldl 
(|l977l ) found fc — 0.05. Assuming a thickness of 
~ Ikpc for the thermal electron s in the disk of 
the g alaxy (the Galactic value, Cordes fc Laziol 
20021) and correcting for the inclination of M33, 
L ~ 1800 pc. Then the extinction corrected Ha 



(EM) map of MSB (jTabatabaei et al.ll2007bh gen- 
erates a distribution of (n,,) across the galaxy with 
a mean value of ~ 0.05 cm^"^ and a most prob- 
able value of ~ 0.03 cm~3 (Fig. [TO]), that is in 
agreement w ith the estimated values in our galaxy 
Cordes fc Laz io 200i) and other near by galaxies 



Krause et aDll989l : lOumke et al.ll2000h . Note that 



a more realistic approach would consider different 
filling factors and electron densities for the thin and 
thick disk of the galaxy. However, because the only 
information we have is a superposition of these com- 
ponents along the line of sight, we are not able to 
distinguish the role of each component. The result- 
ing (ue) and i?tur obtained in Sect l4.3l (Fig. [5]) en- 
able us to estimate DF^ at 3.6 and 20 cm. The left- 
bottom panel in Fig. [Tl] shows the ratio of DP^ at 
20 and 3.6 cm. 

The other Faraday depolarization effect that 
is strong in M33, di fferen tial Faraday rotation , is 
given bv lBurnI (|l966h and lSokoloff et al.l (|l998h as. 



DP„ = sine (2RM, A^), 



(11) 



where for simplicity we assume that the disk of M33 
can be represented as a uniform slab. Using the RM^ 
map in Fig. [31 we estimated DP„ between 20 and 
3.6cm across the galaxy (Fig. [TTl bottom-right). 
As small variations in RM; produce large changes 
in the sine function in Eq. (fTTj) . the resulting DP„ is 
not smoothly distributed among neighboring pixels. 



12 
10 

ST 

— 8 
o 

I 6 
t 4 
2 




mean = 0.046 
stddev = 0.029 



Thermal electron density 

Fig. 10. Histogram of the thermal electron den- 
sity (rie) (cm~^) distribution across M33, derived 
from a n extinction corrected Ha (Tabatab aei et ahl 
l2007b ). The mean and standard deviation (stddev) 
of the distribution are given also. 



Qualitatively, both kinds of Faraday depolar- 
ization contribute to the observed depolarization 
in M33. The global phenomenon, the north-south 
asymmetry, is visible in both DP„ as weaker de- 
polarization (DP20/3 ~ 1) in the north and also 
in DP J. as stronger depolarization (DP20/3 ^ 0) 
in the south part of the central region. However, 
locally e.g. at the positions of HII complexes and 
the southern spiral arms, DP^ could explain the 
observed depolarization. The contributions of DP„ 
and DP,, vary region by region. A more quantita- 
tive comparison requires a combination of DP„ and 
DPr across the galaxy, but this needs a detailed 
modelling of depolarization along with distribution 
of the filling factors / and fc, the pathlength L, and 
the turbulent scale d across the galaxy. 

In the south of M33, a strong turbu- 
lent condition was already indicated from the 
HI line- widths being large r than in the north 
( Deul fc van der HulstI 1987 ). which could be con- 
nected to the high starformation activities in the 
southern arms particula rly in the main arm I S 
( Tabatabaei et al.ll2007ct ). Hence, we conclude that 
the highly turbulent southern M33 along with a 
magneto-ionic medium containing vertical regular 
magnetic fields, reduce the degree of polarization 
of the integral emission from the southern half and 
cause the wavelength-dependent north-south asym- 
metry in polarization (or depolarization). 

6. Discussion 

6.1. Vertical magnetic fields 

The model regular magnetic field described in 
Table [2] and shown in Fig. [7] has a vertical compo- 
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Fig. 11. Top: observed depolarization DP20/3 between 3.6 and 20 cm (Eq. [9]). Bottom-left: estimated 
depolarization DP20/3 due to dispersion in Faraday rotation (Eg. llOp . and bottom-right: estimated depo- 
larization DP20 /3 due to differential Faraday rotation (Eq. [TT|) . The straight line shows the minor axis of 
M33. 
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nent in each ring. In the inner ring, the combination 
of the modes m — zO and stronger m — zl produces 
a sinusoidal vertical field that is strongest near the 
major axis: pointing away from us at ~ 30° and 
towards us at 6* ~ 210°. In the outer ring the ver- 
tical field is also strongest near the major axis, at 
e ~ 10° (directed away) and 6 ~ 190° (directed 
towards). 

The presence of a vertical component to the 
regular magnetic field was already indicated lo- 
cally by our rotation measure maps (Sect. [3|). 
However, the large scale vertical field required by 
our fits requires a global origin. The strong 
along the major axis together with the large line 
of sights through the magneto-ionic medium on 
the eastern and western minor axis (Sect. [3]) sug- 
gest that warp may play a role. In other words, 
the 'vertical' field that we i dentify may be due 
to t he seve re war p in M S B (iRogstad et al. 1976 : 
Reakes fc Newtoiil 19781 Sandage fc Humphrevd 



19801: ICorbein &: Schneideil Il997l) The inner HI 
disk investigated bv^Rog stad etaA. i 1976h . shows a 
warp beginning at a radius of ~ 5 kpc with a change 
in the inclination angle of 40° at 8 kpc. The warp of 
the optical plane begins as close to the center as the 
first arm system at 2 kpc (the center of our inner 
ring) , with a change in the arm inclination of > 15° 
at 3 kpc and 25° at 5 kpc (ISand age &: Humphreys! 
1980>) . The model in Sect. 14.21 assumes a constant 
inclination of i = 56° but in a strongly warped disc i 
varies with radius and azimuth. In this case, even if 
Bj-cg only has compo nents in the warped disk plane 
Bd, as for e.g. M51 (iBerkhuiisen et al.[ll997l), M3 1 



( Fletcher et all 120041) . and NGC6946 llBeckl 120071) 
there will be an apparent vertical component Bz 
(as well as an apparent disk parallel component B^) 
with respect to the average disk-plane. 

The ratio B^/B^ — tanz^, where iw is the 
warp inclination. So for Wi ~ 15° in the inner 
ring Bz/Bd ~ 0.3 and in the outer ring Wi ~ 25° 
gives Bz/Bd — 0.5. Our model field in Table [5] has 
Bz/Bd < 2 ± 1 in the ring l-3kpc and Bz/B^ < 
1.0 ± 0.4 in the ring 3-5kpc. This indicates that, 
in the outer ring, the vertical field could be mainly 
due to the warp. However, a real vertical field of a 
broadly comparable strength to the disk field can 
exist in the inner ring. 

6.2. Magnetic and spiral-arm pitch angles 

The pitch angles of the horizontal component of 
the regular magnetic field are high: pB = 48° in 
the ring 1-3 kpc and ps = 42° in the ring 3-5 kpc. 
These magnetic field pitch angles are however lower 
than the pitch angles of the optical arm segments 



identified bv lSandage fc Humphrey^ (|l980l ). which 
are typically p^ = 60°-70°. A combination of shear 
from the differential rotation, producing an az- 
imuthal magnetic field with pB — 0° , and compres- 
sion in spiral arm segments, amplifying the compo- 
nent of the field parallel to the arms pb — 65°, may 
be responsible for the observed pb — 40°. However, 
this type of alternate stretching and squeezing of 
the field could not produce the m — azimuthal 
mode that is found in both rings, unless the pre- 
galactic field was of this configuration. The pres- 
ence of a significant m = azimuthal mode of -Brcg 
can be explained if a large-scale galactic dynamo 
is operating in M33: the axisymmetric mode has 
the fastest growt h rate in disk dynamo models (e.g. 



Beck et al.l 119961 ) . This does not mean that a dy- 



namo is the origin of all of the regular magnetic 
field structure in M33. In particular it would be 
a strange coincidence if the large pb, higher than 
the typical pb in other disc galaxies by a factor of 
^ 2, is not connected to the open spiral arms with 

Pa ^ 65°. 

A rough estimate of the magnetic field pitch an- 
gles expected due to a simple mean-field dynamo 
can be obtained by considering the ratio of the 
alpha-effect — parameterizing cyclonic turbulence 
generating radial field Br from azimuthal Bo — 
to the omega-effect — describing differential rota- 
tion shearing radial field into azimuthal. This can 
be written as ( Shukurovll2004j ) 



tan Pb = 




(12) 



where a is a typical velocity of the helical turbu- 
lence, h is the scale height of the dynamo active 
layer and G — R dO, / di? gives the shear rate due to 
the angula r velocity VI. Using t he HI rotation curve 
derived bv lCorbelli fc Saluccil (l2007l) a ^ Ikms"^ 
as a typical value (jRuzmaikin et al.lll988r ). and an 
HI scale height of 250pc at i?=2kpc increa sing 
steadily to 650pc at i?=5kpc ( BaldwinI 1981 ) we 
obtain approximate pitch angles of pb — 20° and 
Pb — 15° for the rings 1-3 and 3-5 kpc, respectively. 
These are only about 1/2 to 1/3 of the fitted pitch 
angles of the m=0 modes. Models specific to M33, 
which allow for dynamo action as well as the large 
scale gas-dynamics of the galaxy, are required to 
understand the origin of the large pe as well as the 
vertical component of the regular magnetic field. 

6.3. Energy densities in the ISM 

The energy densities of the equipartition magnetic 
fields in the disk {B^/8it and BI^/Stt for the to- 
tal and regular magnetic fields, respectively) are 
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Fig. 12. Energy densities and their variations with 
galactocentric radius in M33. 



shown in Fig.[T21 The thermal energy density of the 
warm ionized gas, |(ne)fcT'ej is estimated from the 
Ha map assuming Tg ~ 10^ K (see also Sect. [SJ. 
Assuming the pressure equilibrium between the 
warm and hot ionized gas with Tg ~ 10 ^ K and 

Ferrier^ 



an el ectron density of ~ 0.01 (rie) (e.i 
20011 ) . the energy density of the hot ionized gas is 
about the same order of magnitude as the warm 
ionized gas energy density. For the neutral gas, we 
derive the energy density of |(n)fcT using the av- 
erage surface density of total (molecular -I- atomic) 
gas given bv ICorbellil ( 20031) and an averag e tem- 
perature of T ~ 50 K (jWilson et al.lll9971) . The 
warm neutral gas with a typical temperature of 
~ 6000 K has roughly the same thermal energy den- 
sity as the cold neutral gas, d ue to a ~ 100 times 
smaller density (e.g. F crric rea 120011 ). Assuming a 
constant scale height of t he disk of lOOpc (as used 
for NGC6946 lBeckll2007t ). we obtain a gas density 
of (n) ~6cm~^ at R=lkpc (that is about 8 times 
smaller than the corresponding value in NGC6946) 
to ~2cm^'^ at R=5kpc (about 3.5 times smaller 
than that in NGC6946jl. The total thermal energy 
density shown in Fig. [12] includes the contribution 
of the warm, hot ionized and the cold, warm neu- 
tral gas. Figure. [12] also shows the kinetic energy 
density of the turbulent motion of the neutral gas 
estimate d using a turbulent veloci ty of 10km s~^ 
(jCorbcl h fc Schne ider 1997; Milosa vlievi^l2004[ ). 

Generally, the energy densities of all compo- 
nents are about the same order of magnitude as in 
the Milky Way (Cox 2005), but o ne order of mag- 
nitude smaller than in NGC6946 (|Beckll2007t) . The 



small thermal energy density compared to the total 
magnetic field energy density shows that the M33's 
ISM is a low f3 plasma (/? is defined as the thermal to 
magnetic energy density ratio, similar results were 
found in th e Milky Way and NGC6946 Cox 200l 
Beckll2007l) . On the other hand, the small thermal 
energy density compared to that of the turbulent 
motions indicates that turbulence in the diffuse ISM 



is supersonic (as in NGC6946 lBeckll2007t ). This is 
in a greemen t with 3-D MUD mo dels for the ISM 
IdejLvillcz fc Breitschwerdtll2007l ). 

The energy densities of the total magnetic field 
and turbulent gas motions are about the same. 
This hints to the generation of interstellar magnetic 
fields from turbulent gas motions for R<8kpc in 
M33. For a comparison at larger radii, detection of 
the magnetic fields is required using deep surveys of 
Faraday rotation of polarized background sources. 
This seems to be promising as the scale l ength of 
the tot al magnetic field (~24kpc, Tabatabaei et al 



2007b^ is much larger than the present detection 
limit of the radio emission. 

The energy density of the regular field is about 
3-7 times smaller than that of the total field. 
Furtheremore, it shows more variations with radius 
than the total field strength, with a maximum in- 
crease at 4.5 < i? < 5.5 kpc. 

6.4. Comparison with M31 

Comparing magnetic fields in our neighbors, M33 
and M31, is instructive specially because a simi- 
lar method was used to find out the 3-D magnetic 
field structure in both galaxies. It is interesting to 
see that the total magnetic field strength is about 
the same in M33 and M31 (~ 7 fiG), but the reg- 
ular field strength in M33 is about half of that in 
M31 (~ AfiG). Furthermore, in contrast to M33, 
M31 has a disk plane parallel regular field without a 
vertical component. In other words, the large-scale 
magnetic field is well-ordered and strong relative to 
the small-scale field in M31. The regular magnetic 
field is fitted by a d ominant m=0 mode in M31 
(jFletcher et al.l 120041 ) that is much stronger than 
that in M33. 

The larger pitch angles of the horizontal mag- 
netic field in M33 than in M31 is not due to a 
smaller shear in M33. Prom their rotation curves 



^ Hence, the radial profile of the energy density of the 
neutral g as is much fl atter in M33 than in NGC6946 (see 
Fig. 5 in lBeckll2007h 



(jCarignan et al.ll2006t[Corbeni fc Saluccill2007^ ■ the 

shear rate at the relevant radii is larger in M33 
(> lOkms-^kpc^i at 1 < i? < 5 kpc) than in M31 
(~6kms-ikpc-i at 8 < i? < 12 kpc). 

The fact that M33 has a higher star formation 
efficiency than its 10-times more massive neighbor, 
M31, may be a clue to the origin of their differences. 
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Strong starformation activities in the inner part of 
M33 could cause vertical distribution of magneto- 
ionic matter and hence the vertical magnetic field. 
Furthermore, stronger turbulence in the interstellar 
medium can be generally caused by high starforma- 
tion rate increasing the dynamo alpha-effect and 
hence providing large pitch angles of the horizontal 
magnetic field. 

7. Summary 

The distributions of linearly polarized intensity and 
polarization angle at 3.6, 6.2, and 20 cm along with 
the maps of nonthermal intensity and nonther- 
mal spectral ind ex (obtained from the n ew sepa- 
ration method, Tabatabaei et al. 2007b[) yielded 
high-resolution distributions of RM, nonthermal de- 
gree of polarization, and Faraday depolarization in 
M33. Furthermore, we derived the 3-D structure of 
the regular magnetic field by fitting the observed 
azimuthal distributions of the polarization angle 
within two rings of 2 kpc width in the radial range 
1 to 5 kpc. The main results and conclusions are as 
follows: 

1. The average nonthermal degree of polarization 
is Pnth —10% (at 3.6 cm) for R < 7.5 kpc and 
> 20% in parts of the magnetic filament. Due 
to Faraday depolarization P„th decreases to 
~ 6% at 20 cm. 

2. The intrinsic Faraday rotation shows larger 
small-scale variations and weaker correlation 
with PI in the south than in the north of 
M33. The higher starformation activity in the 
southern arms could increase the turbulent 
velocities of interstellar clouds and disturb the 
regular field configuration. On the other hand, 
a good correlation between RM^ and PI in the 
magnetic filament in the north-west of M33 
shows that here the magnetic field is mainly 
regular. 

3. The average equipartition strengths of the total 
and regular magnetic fields are i?tot — 6.4 jiG 
and Srcg — 2.5 /^G for R < 7.5 kpc. The regular 
magnetic field strength is higher within the 
ring at 4.5 < R < 5.5 kpc, which contains 
the magnetic filament that has a maximum 
regular field of i3iog — 6.6 /iG. Strong turbulent 
magnetic fields ( Btm > 7 ^G) occur in the 
extended central region and the arms I S and 
part of II S. 

4. The 3-D structure of the regular magnetic 
field can be explained by a combination of 



azimuthal modes of m = + zO + zl in the ring 
1-3 kpc and m = 0-|-l + zl in the ring 3-5 
kpc. The horizontal magnetic field component 
follows an arm-like pattern with pitch angles 
smaller than those of the optical arm segments, 
indicating that large-scale gas-dynamical effects 
such as compression and shear are not solely 
responsible for the spiral magnetic lines. The 
significant axisymmetric mode (m=0) in both 
rings indicates that galactic dynamo action is 
present in M33. 

5. The presence of vertical magnetic fields, shown 
by the best-fit model of the 3-D field structure 
(zl) and indicated by the Faraday rotation 
distribution across the galaxy, is possibly due to 
both global (e.g. M33's warp or interaction with 
M31) and local (e.g. starformation activities, 
Parker loops) phenomena. The warp can better 
explain the origin of the vertical field in the 
outer ring (3-5 kpc). 

6. In the southern half of M33, an excess of dif- 
ferential Faraday rotation together with strong 
Faraday dispersion seem to be responsible for 
the north-south asymmetry in the observed 
depolarization (which is wavelength dependent). 

7. The energy densities of the magnetic field and 
turbulence are about the same, confirming the 
theory of generation of interstellar magnetic 
fields from turbulent gas motions. Furthermore, 
it seems that the ISM in M33 can be charac- 
terized by a low plasma and dominated by 
a supersonic turbulence, as the energy densities 
of the magnetic field and turbulence are both 
higher than the thermal energy density. 
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